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We study the one-dimensional quantum critical spin systems with the sine-square deformation, 
in which the energy scale in the Hamiltonian at the position x is modified by the function fx = 
sin^ [x'(^ ~ l)]' "*^here L is the length of the system. By investigating the entanglement entropy, 
spin correlation functions, and wave-function overlap, we show that the sine-square deformation 
changes the topology of the geometrical connection of the ground state drastically; Although the 
system apparently has open edges, the sine-square deformation links those ends and realizes the 
periodic ground state at the level of wave function. Our results propose a new method to control 
the topology of quantum states by energy-scale deformation. 
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Introduction: Topology is one of the most fundamen- 
tal concept in physics. It rules the connectivity of local el- 
ements of the system and governs how physical objects - 
particles, excitations, and information - propagate. Nor- 
mally, the topology of a system is fixed once the spatial 
geometry of elements is given. Search for other paths 
to the control of topology of the system is a challenging 
problem. 

In a finite system, a boundary condition determines 
the topology of the geometrical connection of quantum 
state and affects crucially the properties of the system. If 
the system has open edges, they usually induce boundary 
oscillations such as Friedel oscillation. While the bound- 
ary oscillation contains important information such as 
the Fermi momentum, it is often regarded as an obstacle 
to mask the bulk properties. One simple way to remove 
it is to employ the periodic boundary condition, how- 
ever, there has also been another attempts, called the 
smooth boundary condition, to suppress the boundary 
effects by turning off the energy scale of local Hamilto- 
nians smoothly around the open edges. [l|, Q The latter 
has proven to be useful when the open system is favored, 
e.g., for an efficiency of numerical methods such as the 
density-matrix renormalization group (DMRG) method. 

Recently, a new scheme of the smooth boundary con- 
dition, which we call the sine-square deformation (SSD), 
has been proposed as an efficient way to suppress the 
finite-size and open-boundary effects. [3| In the system 
with SSD, the energy scale in the Hamiltonians are mod- 
ified according to the function. 
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where x is the position of the local term and L is the 
length of the system. In Ref. H, Gendiar et al. applied 
the SSD to the one-dimensional (ID) free fermion sys- 
tem with open boundaries. They then showed that the 
SSD removed boundary effects successfully and resulted 
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FIG. 1: (a) Rescaling function of the SSD. (b) Bond 
strength (S^Sj+i) for the XXZ chain © with L = 80 and 
(A, A'f) = (0.5,0). Squares and circles represent the data for 
the chain with SSD and the uniform open chain, respectively. 



in position-independent one-point functions such as the 
bond strength and particle density in the ground state. 
Since the spatial profiles of these quantities were nearly 
completely flat, the observation raised a natural question 
of what happened in the ground state of the system with 
SSD. This is indeed the motivation of the present study. 

In this paper, we study the SSD in several ID quan- 
tum spin systems. Using the DMRG and exact diagonal- 
ization methods, we study numerically the entanglement 
entropy (EE), correlation functions, and wave- function 
overlap in the systems with SSD. We then show that the 
ground state of a critical system with SSD is equivalent 
to that of the uniform periodic system; the SSD changes 
the topology of the critical ground state drastically, from 
an open chain to a periodic ring. The result opens the 
possibility to control the topology of quantum states by 
the energy-scale deformation even in the case that the 
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geometrical shape of the system is fixed. 

Sine-square deformation: The SSD introduces a spa- 
tial modulation of energy scale by applying the rescaling 
factor fx [Eq.(IT|)] to the local Hamiltonian at the position 
X. For example, the model Hamiltonian of the spin-1/2 
antiferromagnetic XX Z chain with SSD is given by 

L-l 
L 

-hY^f^^h (2) 

where we have introduced the magnetic field h which 
induces magnetization M per spinjj]. Hereafter, we con- 
sider the case of even L otherwise mentioned. The energy 
scale of the local Hamiltonians thus decreases smoothly 
as being closer to the boundaries and eventually vanishes 
at the open ends, as shown in Fig. [U (a).0 

Figure[T] (b) gives the DMRG data of the bond strength 
(S'JS'j^+i) in the ground state of the XX Z chain ^ with 
and without SSD. The data clearly show that the SSD 
eliminates the Friedel oscillation seen in the uniform open 
chain almost completely Q. We will demonstrate below 
that the SSD is not only an efficient measure to sup- 
press the boundary effects but also a device to drastically 
change the topology of the ground-state wave function. 

Entanglement entropy: We first investigate EE in the 
ground state of the ID systems with SSD. We consider 
EE for a subsystem O of the left I spins, 

5(0 = -Tro[p(01np(/)], (3) 

where p{l) is the reduced density matrix for VL. For the 
ID critical uniform systems, EE is known to take a uni- 
versal form, 0-11] 

5(0 = sln[.g(0] + const., (4) 

where g[l) = ^sin(^). The slope s is determined by 
the boundary condition; s = c/3 for the periodic system 
while s = c/6 for the system with open boundaries, where 
c is the central charge. Namely, the slope s divided by 
c/6 gives the number of "cuts" of the ID critical state 
between the subsystem H and the environment il. 

Figures [H (a) and (b) show the DMRG data of EE in 
the XX Z chain with SSD [Eq. (0)]. EE in the uniform 
open chains is also shown for a comparison. Remarkably, 
EE in systems with SSD has the slope s = c/3, which is 
twice as large as that in the uniform open systems. This 
means that the ground state of the system with SSD has 
two cuts between the left and right subsystems, and 
although the lattice has seemingly only one cut. In ad- 
dition, the boundary oscillation, which is pronounced in 
the uniform systems, is removed by the SSD. The results 
suggest that, although the system apparently possesses 
the open edges, the SSD connects the open ends of the 
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FIG. 2: Entanglement entropy S{1) as a function of g{l) = 
(L/7r)sin(7rZ/L) for i = 80 and (a) XXZ chain with 
(A,M) = (0.5,0), (b) XXZ chain with (A,M) = (1.0,0.25), 
(c) J1-J2 chain with (J2/Ji,M) = (0.5,0.125), (d) J1-J2 
chain with {J2/ Ji, M) = (0.6,0.4), and (e) two-leg ladder 
with (Jx/J||,M) = (1.0,0.25). The central charge c = 1 for 
the models in (a)-(c) and (e), while c = 2 for (d). Squares 
and circles represent the data for the open system with SSD 
and the uniform open system, respectively. Solid and dotted 
lines show the slope of c/3 and c/6, respectively, (f) Shape of 
the subsystem f2 for which S^V) is calculated. 



ground state effectively and the state becomes periodic^ 
having two cuts between f2 and fi. 

We have also examined EE in the other models, the 
antiferromagnetic Ji- J2 chain and two-leg ladder systems 
under magnetic field. The Hamiltonians are given by 

L-l L 

^Jl-J2 = X! X! Jnfj+^Sj ■ Sj+n - h'YfjSj, (5) 
j = l n=1.2 j=l 
L-l 

j=l n=l,2 

L L 

In Figs. [2] (c)-(e), we present the DMRG results of EE 
for the subsystem f2 of left I sites/rungs [sec Fig. [5] (f)]. 
The models in Figs, ^c) and (e) arc in critical phases 
with c = 1, while the model in Fig.HId) has c = 2[i-[l3]. 
It is found again that the slope of EE is doubled by the 



3 




A 



V 



0.4 



0.2 



XXZ chain, A=1.0, M=0 _ 



XXXXXXXXXXXXXXXXXXXXXXXI 



xxxxxxxxxxxxxxxxxxxxxxl 

XXXXXXXXXXXXXXXXXXXXXHI 

XXXXXXXXXXXXXXXXXXXXBHI 
XXXXXXXXXXXXXXXXXXHHIHI 

xxxxxxxxxxxx 



r=l 



r=2 
r=3 
1=4 
r=6 
r=12 



10 

i '=j+r-L j 



20 

j '=]+>■ 



(b) 



FIG. 4: (a) Spin correlation function (— 1)' (Sj ■ Sj/) with 
J = j+r (mod L) in the XXZ chain for L = 24 and (A, M) = 
(1.0,0) as a function of j and r. Symbols show the data for 
the open chain with SSD; crosses represent the correlations 
between the sites j and j' = j + r (pairs "within" the chain) , 
while squares are those between j and j' = j + r — L (pairs 
"across" the edges) . Lines show the values of the correlations 
in the uniform periodic chain, (b) Schematic picture to show 
the two sites {j,j') at the "distance" r. 



FIG. 3; (a) Spin correlation functions {S^S",) (a = x,z) in 
the XXZ chain for L = 80 and (A, M) = (0.5, 0) as a function 
of the distance \j — where the sites are selected as 

j = L/2-[r/2] andj' = L/2+ [(r + 1)/2]. Squares and circles 
represent the DMRG data for the open chain with SSD and 
the uniform open chain, respectively, while the lines show the 
analytic result for the uniform periodic chain, (b) Schematic 
picture to show the relation between the pairs {j, j') in the 
open chain with SSD and in the periodic chain. 



SSD. We note that, for the J1-J2 chain with large J2/J1, 
sizable boundary osciUations is observed in EE and one- 
point functions (not shown) even in the system with SSD. 
This is presumably attributed to an effective boundary 
field which can not be eliminated completely by the SSD 
in Eq. (in])-^143 However, we emphasize that even in that 
case the doubled slope of EE is observed, which suggests 
that the SSD works also for those models to lead to the 
topology change of the ground state. We thus conclude 
that the change of slope of EE is not peculiar to a specific 
model but a general outcome of SSD when applied to a 
critical model. 

Correlation functions: We next investigate the two- 
spin correlation functions. Here, we consider the spin-1/2 
XXZ chain in the critical regime, for which the asymp- 
totic forms of the correlation functions are known to be 
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where Q = 2ttM. The exponent 77 and the amplitudes 
Aq, Af, and Af have been obtained as a function of A 



and M.[lll[l5l-[l7t Figure [3] shows the DMRG results of 
the ground-state correlation functions in the XXZ chain 
@ with SSD. We also plot the DMRG data for the uni- 
form open chain as well as the analytic result for the 
uniform periodic chain, the latter is obtained by replac- 
ing r in Eqs. ([7]) and ^ with ^ sin I ^ . As seen in 
the figure, the results for the open chain with SSD agree 
almost completely with those for the periodic chain. 

Figure [4] (a) shows the ground-state correlation func- 
tion {Sj ■ Sj') in a small system calculated by the exact 
diagonalization. The data are plotted as a function of the 
position j and the "distance" r = min(|j — L— |j— j'|) 
[see Fig. |4] (b)]. We observe again that the correlations 
in the open chain with SSD are in excellent agreement 
with those in the uniform periodic chain; The results are 
independent of the position j and, more remarkably, the 
correlations between the sites j and j' = j + r — L, which 
locate at the distance r across the open ends, have the 
same value as those in the periodic chain, [l^l We have 
observed the same phenomena as those in Figs. |3]and[3] 
for several parameter sets of (A,il/). The results indi- 
cate that correlation functions, and presumably all ob- 
servables, in the ground state of the systems with SSD 
become equal to those in the uniform periodic systems. 

We note that for the two-leg ladder with zero magneti- 
zation M = 0, which has an energy gap above the singlet 
ground state, the spin correlation decays exponentially 
even in the systems with SSD and no recovery of the cor- 
relation between edge spins is observed. This suggests 
that the SSD does not work for the spin-gapped systems. 

Wave Functions: Finally, we discuss the overlap of 
the ground-state wave functions. Using the exact di- 
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agonalization method, we have calculated the ground- 
state wave function |vssd) in the XX Z chain ([2]) with 
SSD for L < 24 and several sets of (A,M), and com- 
pared it with the ground-state wave function |vpBc) of 
the uniform periodic chain. We have then found that 
the overlap of those ground-state wave functions is very 
close to unity; the deviation from the unity is at most 
|1 — (vssd|vpbc)| ^ 10^"^ and exactly zero within the 
numerical accuracy of 10^^'* for the XX case (A = 0). 
The result indicates that the ground states |vssd) and 
|vpBc) are equivalent at the level of wave function. [ist 

We note that the equivalence of the ground-state wave 
functions is not trivial even in the case of the XX chain 
[Eq. ^ with A = 0] . Through the Jordan- Wigner trans- 
formation, the XX chain is mapped into the free fermion 
system and the one-particle eigenstates of the periodic 
chain is the simple plane waves. On the other hand, the 
Hamiltonian of the open chain with SSD is not trans- 
lationally invariant and its one-particle eigenstates are 
distinct from the plane waves. Nevertheless, when and 
only when the fermions are filled up to the Fermi level, 
the Slater determinants of the two ground states become 
equivalent. This means that the excitation spectrum and 
dynamics of the system with SSD are in general different 
from those in the periodic system. 

Concluding remarks: In summary, we have studied 
the SSD applied to ID critical spin systems. From nu- 
merical analyses of the entanglement entropy, correlation 
functions, and wave-function overlap, we have shown that 
the ground state of the open system with SSD is equiva- 
lent to the one of the uniform periodic system. 

We note that our finding that the SSD realizes the pe- 
riodic ground state is not restricted to a specific model 
but a generic feature of SSD. We have found the change 
in the slope of the entanglement entropy for several spin 
systems and the suppression of boundary effects by SSD 
has also been observed in the free and interacting fermion 
systems Only a condition required is that the sys- 

tem should be critical. This may suggest that the result 
can be understood in a theory applicable to general crit- 
ical systems, such as the conformal- field theory. Investi- 
gating the effects of SSD on the low-energy excited states 
would reveal clues about the mechanism of SSD. 

The results of the present study offer a new scheme 
to modify and further control the topology of the quan- 
tum states by the energy-scale deformation, even under 
the condition that the spatial shape of the system is 
fixed. The approach might be applicable to real systems 
such as ultracold atoms, for which the spatial modulation 
of interatomic interactions has been demonstrated. [20l| 
Searching for an energy-scale deformation to yield the 
other topological change of the ground state must be an- 
other interesting problem. 

We thank Shunsuke Furukawa, Akira Furusaki, Andrej 
Gcndiar, Hosho Katsura, and Tsutomu Momoi for stim- 
ulating discussions. T.H. was supported in part by the 



MEXT and JSPS, Japan through Grand Nos. 21740277 
and 22014016. 



[1] M. Vekic and S. R. White, Phys. Rev. Lett. 71, 4283 
(1993). 

[2] M. Vekic and S. R. White, Phys. Rev. B 53, 14552 (1996). 
[3] A. Gendiar, R. Krcmar, and T. Nishino, Prog. Theor. 

Pliys. 122, 953 (2009); Prog. Theor. Phys. 123, 393 

(2010). 

[4] Unhke the usual uniform systems, the Zeeman term in 
the Hamiltonian with SSD does not commute with the 
exchange term and the relation between the field h and 
the ground-state magnetization M is not trivial. In our 
calculation, we use the value of h{M) obtained from the 
Bethe ansatz or the DMRG analysis for the correspond- 
ing uniform systems as an approximation of the field to 
realize the ground state with M in the system with SSD. 

[5] To be precise, applying SSD to a periodic chain with only 
nearest-neighbor interactions leads to the same Hamilto- 
nian ((21) since fj^_^_i = 0. In this paper, we call the model 
^ the "open" chain with SSD as it apparently has the 
open ends. 

[6] C. Holzhey, F. Larsen, and F. Wilczek, Nucl. Phys. 

B424, 443 (1994). 
[7] G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. 

Rev. Lett. 90, 227902 (2003). 
[8] P. Calabrese and J. Cardy, J. Stat. Mech.: Theory Exp. 

(2004) P06002. 
[9] R. Chitra and T. Giamarchi, Phys. Rev. B 55, 5816 

(1997). 

[10] A. Furusaki and S.-C. Zhang, Phys. Rev. B 60, 1175 
(1999). 

[11] T. Hikihara and A. Furusaki, Phys. Rev. B 63, 134438 
(2001). 

[12] K. Okunishi and T. Tonegawa, J. Phys. Soc. Jpn. 72, 479 
(2003). 

[13] T. Hikihara, T. Momoi, A. Furusaki, and H. Kawamura, 
Phys. Rev. B 81, 224433 (2010). 

[14] The reason why the SSD fails to eliminate boundary os- 
cillations in the J1-J2 chain can be understood if one 
consider the system with Ji = 0, which consists of two 
Heisenberg chains placed with a shift of a half lattice 
spacing. In the case, we should apply the SSD to each 
chain as in Eq. ([2]). This suggests that for finite J2/J1 
one must modify the form of SSD around the open edges 
from Eq. To find the correct form is left for future 
studies. 

[15] N. M. Bogoliubov, A. G. Izergin and V. E. Korepin, Nucl. 

Phys. B 275, 687 (1986). 
[16] S. Lukyanov and V. Terras, Nucl. Phys. B 654, 323 

(2003) . 

[17] T. Hikihara and A. Furusaki, Phys. Rev. B 69, 064427 

(2004) . 

[18] For odd L, we have calculated the spin correlation func- 
tions and wave-function overlap and found that the 
ground state of the open XX Z chain with SSD is equiv- 
alent to the one of the XXZ chain under the antiperi- 
odic boundary condition; The edge spins are connected 
via an exchange coupling J [- {S^Sf + SlSf) + ASlSi] 
with the sign change of the transverse part, which picks 



5 



up the TT-phase of the antiferromagnetic interaction. [20] R. Yamazaki, S. Taie, S. Sugawa, and Y. Takahashi, 

[19] A. Gendiar, M . Daniska, Y. Lee, and T. Nishino, Phys. Rev. Lett. 105, 050405 (2010). 

larXiv:1012. 14721 



